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Abstract We report a novel method for detecting cloud particles in the atmosphere. Solar radiation 
backscattered from clouds is studied with both satellite data and a radiative transfer model. A distinct feature 
is found in the angle of linear polarization of solar radiation that is backscattered from clouds. The dominant 
backscattered electric field from the clear-sky Earth-atmosphere system is nearly parallel to the Earth surface. 
However, when clouds are present, this electric field can rotate significantly away from the parallel direction. 
Model results demonstrate that this polarization feature can be used to detect super-thin cirrus clouds having 
an optical depth of only -0.06 and super-thin liquid water clouds having an optical depth of only -0.01. Such 
clouds are too thin to be sensed using any current passive satellite instruments. 


Super-thin cirrus clouds with optical depths smaller than -0.3 exist globally [McFarquhar et at., 2000; Sun et at., 
2011b]. These clouds are important to the radiation energy balance of the Earth [Sasse/i and Benson, 2001 ; Sun 
et at., 201 1 b]. They also can affect the remote sensing of aerosols [Sun et at., 201 1 b], surface temperature [Sun 
et at., 201 1 a], and atmospheric gases [Christi and Stephens, 2004], For example, the aerosol optical depth (AOD) 
from NASA's Moderate Resolution Imaging Spectroradiometer (MODIS) data could be overestimated by -100% 
when these clouds exist [Sun et at., 201 1 b]. Failing to detect these clouds, the sea surface temperature (SST) 
retrieved from NASA's Atmospheric Infrared Sounder (AIRS) satellite data is -5 K lower at tropical and midlati- 
tude regions, where occurrence frequency of these clouds is high [Sun et at., 2011a]. Climate models must in- 
corporate these clouds correctly in order to account for the Earth's radiation energy budget. To intercalibrate 
other satellite sensors with NASA's future Climate Absolute Radiance and Refractivity Observatory (CLARREO) 
measurements [Wieticki et at., 201 3], knowledge on these clouds is also important for correcting those sensors' 
measurement errors due to light's polarization [Sun and Lukashin, 2013]. 

Due to uncertainties in surface reflectance, the transparent super-thin clouds generally cannot be detected by 
satellite imagers, like the MODIS and the Advanced Very High Resolution Radiometer that only measure the total 
radiance of the reflected solar light [Ackerman et at., 2008], The resulting data products of many satellite and 
ground measurements are biased by these undetected clouds [Sun eta!., 2011a, 201 1b; Omar et at., 2013]. Using a 
strong water vapor absorption channel such as the 1 .38 pm radiance to exclude the surface and low-layer effects 
can be effective on high cirrus [Gao and Kaufman, 1 995] but may encounter difficulties for atmospheres with low 
water vapor [Ackerman etal., 1998]. The reliability of this method is also questionable if the clouds' optical depth 
is smaller than -0.5, when their backscattered intensity is low [Roskovensky and Ltou, 2003]. In addition, super-thin 
clouds may also exist in the lower layers of the atmosphere where there is ample water vapor. The sensitivity of 
the 1.38 pm channel is weak in this region, hampering detection capabilities. NOAA's polar orbiting High 
Resolution Infrared Radiation Sounder multispectral infrared data are usually used with the C0 2 -slicing method 
for detecting thin cirrus clouds [Wylie et at., 1995; Wylie and Menzet, 1999]. However, for super-thin clouds, this 
requires the radiance of their background atmosphere and surface to be very close to that of the reference clear- 
sky environment, which can be difficult as the terrestrial background changes on spatial and temporal scales. In 
addition, this method is problematic when the difference between clear-sky and cloudy radiance for a spectral 
band is less than the instrument noise [Wylie and Menzet, 1999], as for super-thin clouds. 

Currently, the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation lidar [Winker et at., 2007] is 
the only instrument in orbit that can detect super-thin clouds. However, long-term global surveys of super- 
thin clouds using space-borne lidar are limited by the large operational cost and narrow field of view of this 
active instrument. In this letter, we report a novel method for using passive polarimetric instruments such as 
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Figure 1 . The geometry of the system showing the scattered light from a surface located in the x-y plane. The Sun is located 
on the principal plane (xoz) and over the negative x axis (i.e., the solar azimuth angle is 180°). The direction of the reflected 
light is specified by the unit vector e r , and 6 and </> denote the viewing zenith angle (VZA) and relative azimuth angle (RAZ), 
respectively. In the local right-handed orthonormal coordinate system e r = egxe^,, where eg lies in the meridian plane of 
the reflected light beam. The AOLP of the reflected beam in the direction of e r is the angle between the local meridian line 
and the electric vector E of the linearly polarized light, counted anticlockwise when viewing in the reverse direction of the 
reflected beam. The electric field E of the linearly polarized part of the reflected light oscillates in the plane of eg and e^. 


the aerosol polarimetry sensor that was a part of NASA's Glory mission [ Mishchenko et al., 2007] to detect 
super-thin clouds. 


A space-borne polarimetric instrument can measure the Stokes parameters I, Q, U, and V of the scattered light 
[Mishchenko et al., 2007]. Since the circularly polarized component of radiance reflected by the ocean-atmosphere 
system is negligible (l/«0), the degree of polarization (DOP) and angle of linear polarization (AOLP) can be defined 
in terms of Stokes parameters as [Sun and Lukashin, 2013] 


DOP 


VQ 2 + u 2 

i 
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and 
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AOLP = - tan^ 1 - 
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where a o = 0° if Q> 0 and U>0; a 0 = 180° if Q> 0 and U< 0; and a o = 90° if Q<0. In a right-handed Cartesian 
coordinate system, as shown in Figure 1, the AOLP of the reflected light is the angle between the local meridian 
line in the direction of eg and the electric vector E of the linearly polarized part of the reflected light, measured 
anticlockwise when viewed in the reverse direction of the reflected radiance [Sun and Lukashin, 201 3]. The elec- 
tric field E of the linearly polarized part of the reflected light oscillates in the plane of e 0 and in Figure 1 . 

Natural solar radiation can be polarized by surface reflections as well as by scattering from atmospheric 
molecules and particles. When sunlight propagates through the clear atmosphere and is scattered back to- 
ward the Sun, the resulting signal is nearly unpolarized when the solar zenith angle (SZA) is not larger than 
-40° [Sun and Lukashin, 201 3]. The residual polarization in this direction is due to asymmetries in the system, 
due for instance to preferentially oriented ocean waves or nonzero angles between incidence and observa- 
tion. By considering a longer solar wavelength, such as 670 nm, the contribution of molecular scattering is 
reduced. Moreover, unlike total radiance /, the DOP and AOLP are insensitive to surface roughness and 
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Figure 2. The AOLP at 670 nm as a function of VZA and RAZ from (a) the mean of 24 day PARASOL measurements for dear- 
sky oceans identified by the PARASOL product, and (b) the ADRTM for clear-sky oceans. The 24 days of PARASOL data are 
taken from the first 2 days of each month across 2006. Data are binned in SZAfrom 27° to 30° and in wind speed from 6 m/s 
to 9 m/s, and averaged in 3° bins in VZA and RAZ. In the modeling, the SZA is 28°, the ocean wind speed is 7 m/s, the sea-salt 
AOD is 0.06, and the U.S. standard atmosphere is used. 


absorption by atmospheric water vapor and other gases [ Sun and Lukashin, 2013], which makes the polari- 
zation measurement robust in different environmental conditions, even when the detected components are 
within the lower layers of the atmosphere. In brief, remote sensing atmospheric particulates using polariza- 
tion measurements in the backscattering region can minimize surface, molecule, and absorbing gas inter- 
ferences, thus increasing the sensitivity to atmospheric particulates, like super-thin clouds. 

Because of to the variations in surface reflections and atmospheric profiles, using total reflection intensity to detect 
super-thin clouds is generally difficult. Also, although thin clouds can cause changes in the degree of polarization 
(DOP) of the reflected light [Sun and Lukashin, 2013], the dynamic range of the DOP change is mostly insufficient 
to unambiguously identify super-thin clouds if the background polarized reflection is uncertain. However, another 
fundamental characteristics of light is its angle of linear polarization (AOLP). Figure 2 shows the AOLPs at a 
wavelength X = 670 nm (Figure 2a) taken from the mean of 24 day measurements for clear-sky oceans from the 
Polarization and Anisotropy of Reflectances for Atmospheric Science coupled with Observations from a Lidar 
(PARASOL) product [Deschamps et al., 1994], and (Figure 2b) calculated using the adding-doubling radiative 
transfer model (ADRTM) [Sun and Lukashin, 201 3] for clear-sky oceans. The 24 days of PARASOL data are taken from 
the first 2 days of each month across 2006 and collected in a SZA bin of 27-30° and a wind speed bin of 6-9 m/s, 
and averaged in 3° bins in viewing zenith angle (VZA) and relative azimuth angle (RAZ). In the modeling, the 
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Figure 3. The AOLP at 670 nm as a function of VZA and RAZ from (a) the ADRTM for thin cirrus (Ci) composed of solid 
hexagonal ice columns, (b) the ADRTM for thin cirrus composed of a mixture of typical midlatitude ice particle shapes, 
and (c) the mean of 24 day PARASOL measurements for ice clouds over both ocean and land identified by the PARASOL 
product. In the modeling, the optical depth (OD) of the cirrus clouds is 0.1 at 670 nm, the wind speed is 7 m/s, the sea-salt 
AOD is 0.06, the SZA is 28°, and the U.S. standard atmosphere is used. 


SZA is 28°, the ocean wind speed is 7 m/s, the sea-salt aerosols are of irregular shapes with an AOD of 0.06, and 
the U.S. standard atmosphere [National Oceanic and Atmospheric Administration et a!., 1976] is used. Note here 
that we only display the AOLPs over the RAZ range of 0° to 1 80°. The AOLPs over the RAZ range of 1 80° to 360° 
are symmetrically supplementary angles of the values over the RAZ range of 0°to 180° [Sun and Lukashin, 2013]. 
We find that the AOLPs from the PARASOL data and the model are similar except that at ~8° off the exact- 
backscatter direction (RAZ = ~1 80° and VZA = -20° and -36°), the PARASOL result shows two distinct features in 
the AOLP of the reflected solar light. At these viewing angles, the dominant backscattered electric field in the 
PARASOL data for clear-sky oceans rotates from parallel to perpendicular to the Earth's surface; whereas the 
modeled AOLPs do not show these features. We performed sensitivity studies with the ADRTM and found that 
ocean surface roughness due to winds and atmospheric gas or water vapor absorptions do not produce these 
features in the model outputs. However, this feature does appear when we incorporate a thin layer of cirrus into 
the atmosphere. The ADRTM simulation results of Figures 3a and 3b demonstrate that the presence of a super- 
thin layer of cirrus clouds can reproduce the AOLP features in Figure 2(a). In the simulation of Figure 3a, we in- 
corporate the solid-column ice crystal shapes for the cirrus cloud, which may not be the ubiquitous ice crystal 
shapes in the atmosphere. Thin cirrus clouds are generally a mixture of various ice crystal shapes such as hex- 
agonal columns, plates, droxtals, bullet rosettes, and aggregates. In the simulation of Figure 3b, we incorporate a 
typical size distribution and particle-shape mixture of midlatitude cirrus clouds [Heymsfield and Platt, 1984]. While not 
as sharp as the results of Figure 3a, we can see in Figure 3b that the distinct polarization-rotation features at the near- 
backscattering angles are still significant for clouds with an optical depth (OD) of only 0.1 . These features remain visible 
even down to ODs of -0.06. Our sensitivity studies on typical tropical cirrus produce similar results. Further evidence of 
the polarization-rotation features of ice clouds can be seen in the PARASOL data. Figure 3c shows the AOLPs from the 
mean of the 24 day PARASOL data for ice clouds over both ocean and land identified by the PARASOL product. These 
data are particularly noisy because of limited sampling of PARASOL at several viewing angles for ice clouds at the 
specific SZA. However, we can still see the polarization rotations above and below the exact-backscatter angle. 
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Figure 4. The AOLP at 670 nm as a function ofVZAand RAZfrom (a) the mean of 24 day PARASOL measurements for water clouds 
over oceans identified by the PARASOL product, (b) the ADRTM for water clouds with an optical depth (OD) of 0.1, and (c) the ADRTM 
for water clouds with an OD of 0.01 at670nm over oceans. In the modeling, the wind speed is 7 m/s, the sea-salt AOD is 0.06, the SZA 
is 28°, and the U.S. standard atmosphere is used. 


One mechanism that can produce such polarization-rotation features is the optical glory. The glory is an angular 
region of high intensity that may extend several degrees from the exact-backscatter direction and is accom- 
panied by a strong p-polarization component a few degrees from the exact-backscatter direction. This promi- 
nent p-polarization feature appears to be due to the internal field transmitted into the particle, and spherical 
droplets in water clouds can produce an especially strong glory. Figure 4a shows the AOLPs of the averaged 24 
day PARASOL data of water clouds over ocean. It is clear in this figure that the glory feature is a p-polarized band 
that extends -20° from the exact-backscatter direction (VZA - 28°, RAZ = 1 80°). Figure 4b shows ADRTM results 
for water clouds with an OD of 0.1 , the same OD as the cirrus clouds in the simulations of Figure 3. Since the 
mean OD (-7.5) of water clouds from PARASOL data is much larger than 0.1, the p-polarized band in Figure 4a 
extends larger than that in Figure 4b. Also, model result in Figure 4c demonstrates that subvisual water clouds 
having an OD of only 0.01 still display a prominent polarization feature. This polarization feature can be 
suppressed when particle absorption increases and morphology becomes more irregular [ Muinonen eta!., 2011; 
Volten et al., 2001]. While cirrus clouds composed solely of irregularly shaped aggregated ice particles may not 
portray this feature [Sun and Lukashin, 2013], thin cirrus clouds generally contain significant quantities of simple 
particle shapes such as hollow or solid columns, plates, droxtals, and bullet rosettes [Heymsfield and Platt, 1 984], 
especially where close to cloud top or in cold and dry regions. Additionally, since the glory prominence 
is dependent on particle morphology and absorption, our ADRTM simulations incorporating a dust cloud 
[Volten eta!., 2001] do not result in a noticeable glory feature in the AOLP. Such results are consistent with Flubble 
Space Telescope images taken during the 2003 Mars opposition that also showed p-polarization features, which 
were suggested to result from thin clouds of nucleating ice crystals. These crystals were seen at the forward edge 
of a prominent dust storm and were visible against a desert background, but the dust storm itself did not display 
such p-polarization features [Shkuratov et al., 2005]. This demonstrates that polarization-based detection of 
super-thin clouds can be used over different terrain types and can differentiate such clouds from dust. 

In summary, the angle of linear polarization (AOLP) of scattered sunlight observed in two distinct angular re- 
gions near the exact-backscatter direction rotates from an angle parallel to the Earth's surface to perpendicular 
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to the Earth's surface when cloud particles are present in the atmosphere. A sensitivity study demonstrates that 
the polarization-angle feature of background aerosols is not prominent enough to overcome the background 
polarization; however, this feature is detectable in water clouds having an optical depth of only -0.01 and in ice 
clouds having an optical depth of only -0.06, when the solar zenith angle (SZA) is not much larger than -40° 
over a strong polarization background such as oceans. When over weak polarization land, the SZA can be even 
larger. Since most super-thin clouds are over the tropical and midlatitude regions [Sun et at., 201 1 a, 201 1 b], this 
limit of SZA will not affect the detection of most super-thin clouds using the polarization-rotation features. 
Coupling this method with the total radiance measurements at a wavelength of 1.38 pm would enhance the 
detection of super-thin clouds and thus can significantly improve retrieval of aerosol, water vapor or other 
gases, and surface temperature from satellite data. 
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